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SHELL  PARAMETER  INFLUENCE 
ON  THE 

FAILURE  MODE  OF  SIDS  -  ON  AIR  BLAST  LOADED  CYLINDERS 


The  failure  of  aerospace  structures  subjected  to  side-on  air  blast 
loading  has  received  intensive  investigation  by  several  different  groups 
in  the  last  few  yeai *s.  The  results  of  these  intensive  investigations 
indicate  that  the  main  structures  of  aerospace  vehicles  (cylinders) 
fail  in  one  of  three  modes.  These  mode  shapes  are  hinge  collapse,  wave 
buckling,  and  combination  patterns.  The  hinge  collapse  mode  is  charac¬ 
terized  by  the  formation  of  a  plastic  hinge  line  perpendicular  to  the 
longitudinal  axis  of  the  cylinder.  The  wave  buckling  mode  is  character¬ 
ized  by-  the  formation  of  ripple  Waves  in  the  circumferential  direction 
and  one-half  wave  in  the  longitudinal  direction.  The  combination  pat¬ 
tern  is  characterized  by  the  simultaneous  presence  of  the  previously 
mentioned  modes. 

The  time  history  of  the  actual  formation  of  the  failure  pattern  of 
these  aerospace  structures  has  not  been  determined,  rather  a  criterion 
has  been  established  that  specifies  if  a  structure  will  fail.  This 
criterion  states  that  when  the  pressure  and  impulse  generated  by  the 
blast  reach  the  proper  magnitudes  at  the  position  of  the  structure, 
the  structure  will  fril.  Those  magnitudes  of  pressure  and  impulse  which 
produce  failure  are  call  si  critical  pressure  and  critical  impulse.  The 
tabulation  of  these  values  is  dependent  upon  the  structure  material, 
type  of  explosive,  and  structure  position  from  the  point  of  detonation. 
Once  having  such  a  set  of  values,  the  construction  of  iso-damage  curves 
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is  possible.  These  iso-damage  curves  are  the  result  of  a  graphical 
representation*  on  a  logarithmic  scale,  of  the  critical  values  of  pres¬ 
sure  and  impulse*  The  shape  of  these  curves  resemble  hyperbolas  and 
define  tha  boundary  between  failure  and  non-failure  of  a  given  structure. 

The  cue  disadvantage  of  the  above  mentioned  curves,  is  that  they 
(curves)  do  not  predict  how  a  structure  will  fail,  or  what  is  the  domi¬ 
nant  shell  parameter  in  determining  the  type  of  failure  mode. 

The  first  section  of  this  paper  deals  with  background  material.  The 
second  section  defines  the  problem,  lists  the  assumptions,  and  gives  the 
scone  and  method  of  approach.  The  third  section  developes  the  theory 
used  in. discussing  the  various  failure  modes.  The  next  section  dis¬ 
cusses  the  effect  of  the  shall  parameters  in  determining  the  failure  mode 
and  the  last  section  presents  the  conclusions. 
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IX,  Background 

One  of  the  major  problems  in  the  design  of  structures  is  that  of 
determining  their  useful  life.  In  determining  the  useful  life  span  of 
a  structure  it  is  necessary  to  define  the  failure  modes  of  the  structure 
and  the  reasons  for  failure.  This  needs  to  be  done  so  that  design 
judgements  can  be  made  as  to  how  to  reinforce  a  particular  shell  geom¬ 
etry  to  prevent  its  failure. 

In  extensive  experimentation  conducted  by  Scbuman,  at  Ballistic 
Research  Laboratories,  Maryland  (Ref  13),  the  conclusion  was  reached 
that  when  an  aerospace  vehicle  is  subjected  to  sido-on  Hast  loading  of 
sufficient  magnitude  to  cause  incipient  failure ,  it  will  fail  in  one  of 
three  ways  -  hinge  collapse,  wave  buckling,  or  a  combination  pattern. 
Since  most  aerospace  vehicles  are  primarily  of  cylindrical  shape, 

Schunan  used  as  his  model  cylindrical  shells  with  rigidly  clasped  ends. 

He  surmised  that  the  clasped  ends  would  closely  approximate  the  bulk¬ 
heads  found  between  actual  sections  of  an  aerospace  vehicle.  As  a 
result  of  clasping  the  ends,  he  assumed  that  no  pre-loads  such  as  bend¬ 
ing,  compression  or  torsion  were  present.  In  order  to  construct  iso- 
damage  curves,  Shuman  chose  as  his  failure  criterion  a  permanent  radial 
deflection  of  10$  of  the  original  diameter.  After  determining  if  the 
shell  met  the  failure  criterion,  he  used  a  set  of  curves  given  in 
Reference  9  to  obtain  tee  critical  values  of  pressure  and  impulse.  These 
values  were  then  ctey-ared  to  tee  values  obtained  from  testrvimenied  shells 
and  tee  agreement  was  found  to  be  within  the  experimental  error.  Having 
constructed  tbs  iso^damsge  curves,  Scnuman  was  teen  able  so  predict  at 
what  combination  of  pressure  and  impulse  a  shell  of  given  dimension  and 
material  would  be  expected  to  fail.  The  curves  constructed  were  only 
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for  three  types  of  aluminum  and  one  r&ld  steel.  These  iso-daaage  curves 
made  no  attest  to  predict  the  type  of  shell  failure  mode, 

J,  £.  Greehspon  of  J.  G.  Engineering  Research  Associates  (Ref  10), 
U3«d  an  energy  approach  to  calculate  the  energy  necesisary  for  a  given 
deflection  shape  to  occur.  He  used  elastic  stress  analysis,  aesbrane 
theory*  aid  Von  Hises  yeild  condition  to  determine  an  expression  for  the 
critical  collapse  blast  pressure  aid  a  formula  derived  by  Tiaoshenko 
(Ref  14)  for  critical  buckling  blast  pressure.  Ey  using  the  sssaller  of 
these  two  values,  Greenspon  was  able  to  predict  whether  the  shell  ecl¬ 
ipsed  or  buckled  uider  side-on  air  blast  loading.  Once  the  failure 
mode  «as  determined,  he  preceded  with  a  plastic  analysis  to  calculate 
the  energy  of  deformation.  Than  assuring  that  the  energy  of  deformation 
equals  the  initial  kinetic  energy  of  the  shell  due  to  the  blast  loading, 
which  is  related  to  the  applied  impulse,  he  was  able  to  deters! ns  the 
critical  isqjulse  required  for  shell  failure.  Relating  this  impulse  to 
critical  pressure  through  relationships  developed  by  H,  J,  Goodman 
(Ref  9),  Greenspon  was  able  to  construct  a  set  of  iso-dasage  curves  for 
collapse  arid  buckling.  However  these  curves  still  did  not  determine  the 
role  that  each  shell  parameter  played  is  determining  the  failure  isode. 

D.  1,  Anderson  and  H.  S.  Idndbarg,  Stanford  Research  Institute 
(Ref  2),  developed  a  dynamic  buckling  theory  for  pulse  loaded  cylindri¬ 
cal  stalls.  In  their  development  they  used  three  mathematical  models  to 
describe  shell  response  to  pressures  varying  from  ideal  impulsive  load¬ 
ing  to  step  fueticn  loading*  They  chose  the  tan  ent  modulus  model  for 
impulsive  loadings  which  produced  plastic  flow  buckling,  “mey  chose 
the  s trail  .-eversal  ssode 1  for  intermediate  duration  loadings,  which  pro¬ 
duced  plastic-elastic  buckling.  They  chose  the  elastic  model  for  step 
pulse  loadings,  which  produced  elastic  buckling,  using  these  models  in 
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connection  with  experimental  data,  they  were  able  to  determine  the 
equations  for  the  f and  lies  of  .so-damage  curves.  The  resultant  equations 
vers  sensitive  to  variations  in  geometry*  Young* s  aod’ilus,  and  shell 
density.  The  end  conditions  on  their  shells  were  assiased  to  be  simply 
supported  and  infinite  in  length.  Their  reasoning  will  be  discussed  on 
page  11. 

John  P.  Anderson  are!  James  H»  Woodward  at  Kaaan  Nuclear,  Colorado 
Springs,  Colorado  (Ref  3).  approached  the  case  of  dynamic  hinge  collapse, 
due  to  realistic  blast  pressures,  using  an  energy  analysis  and  the 
assumptions  of  finite  length  and  non- symmetric  side-on  air  blast  load¬ 
ing.  The  particular  energy  method  used  was  the  principle  of  minimum 
potential  energy.  The  pressure  distribution  that  they  used  was  an 
exponentially  decaying  cosine  distribution.  The  equations  resulting 
from  this  approach  are  then  solved  for  a  steeply  supported  cylinder  using 
a  Fourier  Sine  series  which  has  time  varying  coefficients  as  the  assumed 
deflection  shape.  This  deflection  shape  corresponds  to  a  simply  supported 
shell,  A  collapse  failure  was  assumed  to  occur  whenever  the  maximum 
moment,  calculated  from  the  deflection  shape,  equaled  or  exceeded 
Brazier’s  collapse  moment  (Ref  6:109), 

In  a  later  work  by  Anderson  and  Woodward  (Fvi  4),  a  static  hinge 
collapse  theory  was  developed.  This  theory  made  use  of  existing  shell 
theory  and  the  assumptions  that  the  cylinder  is  elastic,  isotropic,  and 
infinite  in  length.  The  final  form  of  their  equations  was  derived  by 
assuming  that  1)  the  deflection  shape  could  be  represented  by  a  Fourier 
Sine  series  with  unknown  coefficients,  2)  the  beam  bending  mode  was 
responsible  for  the  primary  response  of  the  shell,  3)  the  breathing  and 
ovalisaiion  modes  are  secondary,  and  4)  the  shell  is  inextensionai  in 
the  circumferential  direction.  This  last  assumption  made  it  necessary 
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to  modify  equilibrium  equations  to  resove  their  sensitivity  to  this 
condition. 

In  spite  of  the  fact  that  all  of  these  researchers  have  investigated 
the  sane  general  problem,  nasely  the  response  of  thin  walled  cylindrical 
shells  to  transverse  pressures ,  no  results  are  presently  available  which 
can  tell  the  engineer  the  answer  to  two  basic  questions.  First,  which 
shell  parameters  play  the  most  important  role  in  defining  the  iso-damage 
curve  of  a  given  shell.  Second,  what  sode  of  failure  will  result  from 
application  of  a  given  pressure  and  isiulse. 


IH.  Problem 


Problem  Statement 

The  problem  is  to  determine  how  the  shell  parameters  influence  the 
failure  modes  of  right  circular  cylinders*  with  clasped  ends*  due  to 
side-on  air  blast  loading. 

The  shell  parameters  are  length  to  diameter  ratio  (L/D),  diameter 
to  thickness  ratio  (D/h),  density  {/*)#  and  modulus  of  elasticity  (E). 

Scope 

The  pressures  mentioned  in  this  paper  are  all  reflected  pressures 
as  opposed  to  incident  pressures. 

The  blast  loading  on  the  shells  are  cosine  squared  and  exponential 
decay  cosine  distributions.  These  loadings  are  on  the  lateral  surface 
facing  the  explosive  charge.  Thus  only  side-on  loading  is  considered. 

The  coordinate  system  that  will  be  used  and  the  failure  patterns 
that  will  be  discussed  are  illustrated  in  figures  (l)  and  (2). 

The  shells  of  interest  are  those  with  L/D  ratios  of  1  to  5  and  D/h 
ratios  of  100  to  ?00.  These  shells  have  been  subjected  to  pressures 
arid  impulses  that  range  from  impulsive  loading  to  step  function  loading. 
The  end  conditions  of  the  shells  are  considered  to  be  clasped.  The 
shells  are  made  of  1040  steel,  6O0I-T6,  505 2-B38,  and  1100-0  aluminum 
and  are  void  of  internal  or  external  reinforcements. 

Assumptions 

In  addition  to  the  assumptions  made  by  the  researchers  listed  in  the 
background  section,  the  following  assumptions  are  made. 

1.  The  shells  are  void  of  ary  pre-stresses. 
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2,  The  relationship  between  the  assent  and  center  line 
curvature  is  valid  fcr  every  cross  section* 

3*  The  body  forces  are  sere, 

4*  The  breathing  node  is  neglected  and  is  independent 
of  collapse  sode- 

5-  The  deflection  in  the  longitudinal  and  tangential 
directions  are  small  compared  to  the  deflections 
in  the  radial  direction. 


Method  of  Approach 

1.  Develops  the  necessary  theory  to  predict  the  criti¬ 
cal  values  of  pressure  as!  impulse. 

2.  Plot  several  iso-dasage  curves  and  analyse  the  gen¬ 
eral  effect  of  L/D  and  D/h  ratios  cn  the  particular 
failure  obtained.  37  this  analysis  the  sost  issor- 
tant  parameter  can  be  determined. 

3.  Prca  the  equations  developed  in  step  (l),  solve  for 
a  pressure  parameter,  f>r  xtoVe  *  and  graph  this  para¬ 
meter  vs  the  most  important  shell  ratio  while  hold¬ 
ing  the  other  ratio  constant.  By  plotting  this 
pressure  parameter  for  each  theory  -  dynamic  and 
static  buckling  and  collapse  -  the  upper  and  loser 
bounds  for  each  type  of  failure  node  can  hopefully 
be  obtained. 

This  method  depends  on  the  developed  theory  for  the  dynamic  collapse 
case  and  the  solution  of  the  subsequent  expression  for  the  critical  pres¬ 
sure.  The  necessary  theory  for  the  cases  of  dynamic  and  static  buckling 
and  static  collapse  have  essentially  been  developed  and  solved  by  previ¬ 
ous  researchers. 


8 


GAW/MC/6S-13 

the  dynamic  equilibrium  equation  in  the  radial  direction  per  unit  axial 
length,  becomes 

-iS  <» 

uhere  N@(©,t)  is  the  membrane  force 0  Q(©,t)  is  the  shear  force  per  unit 
length,  arai  is  the  side-on  pressure,  positive  as  shown  in  Figure  3. 

By  neglecting  rotary  inertia  and  suimaing  moments  on  the  element  shown  in 


Figure  3* 


I  3M 
Q  '  3© 


m  =M(e,t) 


hl^L 

1  -  vz 


vhare  £  is  toe  tangent  modulus,  I  is  the  moment  of  inertia  of  the 
shell,  aid  2/  is  Poisson* s  ratio.  Defining 


T  = 


q*  = 


p(e;r)  =  p(e,t) 


a(s  -vL) 

Eh 
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Substituting  equations  (2),  (4)  to  (11)  by  use  of  equation  (1)  into 
equation  (3)  yields 


y  d  E-j-  3  y 

~~3  "h  — - - f  ■+■ 

©t2  e  ae 


r  a%  Q^Q -v2)' 
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02) 


dine©  the  shell  is  assumed  inertensional ,  the  effects  of  the  flexual 
modes  on  the  average  circumferential  strain  is  restricted  to  second 
order  terms  (Ref  8:269).  The  assured  radial  def creation  and  pressure 
distribution  are  given  by 

CO 

y 


y(8,T)  =  ^(T)  + 


COS  h© 


03) 


co 

^ - i 


Pf (e,T)  -  P«(T)  +  ^  p'n(T)  COS  h© 


(iH) 


n=i 


Substituting  these  relationships  into  equation  (12),  equating  coeffic¬ 
ients,  and  assuming  W0Cr)  is  negligible  compared  to  unity  (Hef  2:205) 
the  resulting  equations  of  motion  are 


a  «o  .  <5(»  “^  )  _  pi 

— .  -p  -j -  _  ^ 

erz  e 


(i  5) 


a2w  2  T  h2azEt 
— f  4-  {nz-u  f 

3  Tc 


Oq{  i-v2)"1 


wn  ~  0&- 


n  =1,2,3--* 
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For  impulsive  loads,  the  pressure  terms  in  equations  (15)  and  (16)  van¬ 
ish,  Trie  infcial  conditions  are 


uh(  0)  =  0 


&±{0)  -  - IlL 

ar  johc 


where 


i{e)  =  ic  + 


Y_,  Tn  COS'( 


n  =  0,5, 2, 3 
cos'(nQ) 


is  the  impulse  applie.-  to  the  shell. 

Equation  (15) »  which  governs  the  motion  of  the  radial  mode,  is 
solved  in  two  parts,  depending  upon  the  value  of  -  If  we  take  Oi 

and  £  as  the  yield,  stress  and  yield  strain  and  Eo  a.  the  stress- 
strain  modulus  beyond  the  yield  point,  then 


rrr"=  >  os^o-v2^ 


°e  1  -  v 


°e  +  T 


«J0  — ,  (l  W0 


Thus  equation  (15)  becomes 


^  4-  w0  --  po 


0<  wQ<  (1— V  ) 


5TUo 

3TS 


ky0  =-(l-  T-](l-^2)ey+Po  ,  « 


which  can  }.*eadily  be  solved. 

Equation  (lo),  which  governs  the  motion  of  the  flax>irsl  modes,  is 
solved  by  assuming  and  E*  constant  (Ref  2:206).  Applying  condi¬ 
tion  (17)  aid  (13)  tc  equation  (16)  the  solution  is 
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,  /-ri  -  SINH  qn{T  -Ty) 

*h'T)  -  ST  qn 

where  Is  the  time  at  which  yielding  first  occurs  and 

<fc(l  ~  ) 


2  =  (n2  - 


*h 


2  2 
-  rra 


£t 


(2H) 


(25) 


The  most  amplified  flexural  mode,  n,  occurs  at  the  naxiaua  value 
of  which  gives  n  as  the  integer  nearest  to  (Ref  2:214) 

jr - - 1  t - - - » 

/j_  i  -vi)  1  _  /_t3iLi_L  (24 

n  =^V  2  l  a2  Et  /  -  h  V  E+ 

Equation  (26)  shows  that  the  most  amplified  flexural  node  is  directly 

dependent  on  a/h  and  \f  ~ t  ♦  From  equriion  (255 »  the  ratio 

-  2 

determines  the  sign  of  thus  determining  Aether  the  motion  is 

hyperbolic  or  oscillatory.  Hence  is  the  material  property  that 

sets  the  magnitude  and  type  of  motion  for  &  shell  under  1  repulsive  load¬ 
ing  (Ref  2:215) « 


Elastic  Ifedel.  The  deformation  of  shells  in  this  region  is  charae 
terised  by  long  duration  pressures  and  buckle  such  that  their  response 
can  be  assumed  as  elastic.  Since  the  buckling  behavior  is  elastic,  a 
reasonable  extension  of  static  theory  is  made  by  the  inclusion  of  iner¬ 
tia  tarns  (Ref  2:198).  To  determine  the  membrane  force,  the  erd  effects 
are  neglected,  but  it  is  required  that  the  flexural  mode  satisfy  the 
boundary  conditions  (Ref  2:193). 

The  governing  equation  of  motion  is  Donnell 5  s  equations  (Ref  5)  *»*ith 
the  addition  of  inertia  terms.  Thus  the  equation  of  equilibria  in  the 
radial  direction  is 
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„  4  0%/ 

DfV  V  + 


■5“  + 


2N_  a“v  n©  a2V 


_xe  _ 

=*  0x50 


02© 


N 


©  ,  S^W 

+  J°h  -P  =  0 


(27) 


where  ,  «x6t  «g  are  the  membrane  forces  positive  in  compression,  and 


D  - 


Eh 


?  !2{l  -  V2  ) 


(28) 


is  the  flexural  rigidity  of  the  shell,  1L  can  be  expressed  as  (2:20?) 


k.  _  th  H>  a2F 

N"  -  "  “  “  +  —~2 


e  (,-v2)  a 


(29) 


where  F  is  the  stress  function  for  the  membrane  forces  due  to  the 


flexural  action  and  WQ  is  the  uniform  radial  deformation.  The  membrane 


forces  and  N  can  now  be  represented  as  (Ref  2:208) 

JL  Xw 


N..  = 


Ac*- 

o  r 


x  a2  a©2 


N. _ =  - 


32F 


x©  a  0x3© 


DO) 

(31) 


Utilizing  the  compatibility  condition 


V4f  =  # 


02W 

0x2 


(32) 


introducing  the  non-dimensional  variables 
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and  expressing  v  and  p  in  series  fora 


yCfe/j  =wc{T)  +  )  wn{T}cosn©^‘n 


p'ge^)  =  Jpo(T)  V] 


equation  (32)  can  be  represented  as 


4(4  *  4]f  --  ^cos^TTl 


d^l  a§z  +  ae^l 


7LIT) 

n=i 


"  (33) 


Thuj' 


-  y*  Bnwn  cos-h©^m^£ 


where  Bn  are  constants.  This  expression  for  F  assuses  simple  supports 
and  that  tbereexisis  only  one  half-wave  buckle  in  the  axial  direction.. 

These  assumptions  are  based  on  experimental  results  of  dynamic  loadings 

(Sef  2:203-9). 

using  equations  (2$)  to  (31)*  (33)  (38).  and  (46),  the  equations 
of  cation  can  be  separated  into 

i#  +  '*  "  *  (Ml) 


a2 

a  ^ 


*  la  (n 


2i’2  _  3T2! 


i  a  i 


I 


(f-V^ 


w  :  *  • 

—  n  wjW., 
o  •  h 


n*i£p-*- 
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with  initial  conditions  as 

du(°) 

wn(°)  -  ~  0  ,  ^3) 

•  Solving  equation  (4l)  for  'the  static  case  and  substitution  of  the 
results  into  equation  (42)  yields 


h 


u 


n 


~  ?.f  2  .  -TT2]2  ,  (i-^)Cn-4/i4r 

d  jfi  T  - 2  ■  + 


(4H) 


i2 1 


2  ir- 
n  +  — T 

r  i 


2  * 
-n-p 

ro 


?roa 
the  s 


equation  (44)  as  the  denominator  approaches  zero  ,  ^  and  thus 
iritical  pressure,  ( pe' }  cr  ,  is  given  ty 

4 


(Pot  = 


n2 


'  2.n2  .  -n-2,  .  (i-^fnyi) 
a  in  +  —, )  .+  — ? — =5-5 


(t5) 


ior 


/  ?-  1 rzf  1 

ih  +  ~^r)  J 

n  =  12,3,4 

2:210). 

By  not  assuming  the  static  case  for  equation  (41),  the  expression 
;q  beeoaes 


W0  =  p0'  +  A.pmT  4.  B  COST  (q4 

am  substituting  this  expression  into  equation  (42)  results  in  Kaihieu 
type  of  solutions  for  vn  (Ref  5:263),  am  if  takes  at  face  value  pre¬ 
dict  instabilities  at  very  lew  values  of  coefficients  A  and  B. 
Instabilities  of  this  type  require  a  constant  energy  input  which,  in 
this  case,  sash  cose  fro-  the  kinetic  and  strain  energy  of  the  ucifora 
hoop  sode.  Huserical  integration  of  equation  (42)  shows  that  it  require 


19 


xAH/HC/6 8-13 


n  large  number  of  oscillations  o?  the  hoop  mode  before  the  flexural 
amplitudes  build  up  to  the  point  wheru  there  is  sufficient  energy  trans¬ 
fer  to  cause  permanent  deformation  (Hof  5:264).  Experiments  show  that 
buckling  takes  place  with  little  or  no  oscillations  ana  is  essentially 
a  single  growth  to  large  deformations,  thus  the  static  case  which 
results  in  hyperbolic  growth  is  a  valid  assumption  {Ref  2:21a) .  Hence 
eouation  (42)  can  be-  rewritten  as 


2  -I-  "*(<»£  -  wo)  «n  «  p' 


where 


.  ,2  2 1  *,2  .  1 
(nOJn)  =  a  ri  +  ~2 


n rV.  (i  -  v2;  {'|r/ 1 ) 


7T2]2 

I2 


and  the  solution  is 


dT 


wn  = 


nVfJ-WK 


*  "5 

jmh  i?s/p;-w|Tj 


for  Po>0Jn 


The  most  amplified  mode  of  equation  (49)  corresponds  to  the  s-axi— 

'If 

^1  of  h^(da-ic)  *  For  T"  ssd  "tno  ssconQ  tsrs  of 

*  v  u  « — 

(48)  small  compared  to  the  first  term,  (Ref  2:219)*  the  sost  amplified 


noce  as 


!  <3  \ 


--V'V^  [Tj 


Equation  (pO)  shows  that  the  most  amplified  liexural  vaas,  n,  is 

.  .  3fz 

directly  proportional  to  \s/h)  -  Thus  this  ratio  seems  to  dominate 

elastic  behavior. 
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Strain  Reversal  Model.  In  (Ref  2i2j3)  it  is  stated: 

pTne  strain-reversal  scdel  to  be  described  in  this  sec¬ 
tion  predicts  threshold  buckling  loads  in  the  portion  of  the 
quasi-isss-ulsive  load  region  in  which  neither  the  tangeni- 
rsodulus  nor  the  elastic  nod  el  is  applicable.  This  is  done  or 
taking  into  account  strains  beyond  the  valid  range  cf  the  elas¬ 
tic  node!  and  allowing  strain  reversal  and  reloading  which  is 
beyond  the  valid  range  cf  the  tange'rbesodulus  sod  el,  which 
treats  only  the  initial  inward  notion  of  the  boon  node.  The 
central  conclusion  fros  the  solutions  using  the  strain  rever¬ 
sal  Hodel  is  that  the  curves  established  by  the  simpler  models 
exhibit  the  essential  characteristics  found  fron  the  sore  cos— 
slex  theory.3 


The  interpretation  cf  the  above  discussion  is  that  a  sisale  exten¬ 
sion  of  the  tangent-modulus  model  and  the  elastic  sodel  into  the  region 
governed  by  the  strain  reversal  esdel  vdl.1  provide  solutions  consistent 
wife  experissnial  results.  This  is  illustrated  on  Figure  ?,  page  33  of 
this  study.  Hence  no  further  discussion  of  this  model  will  fee  presented 


is  this  study. 


Static  Budding  Theory 


In  this  melding  analysis,  the  stability  criterion  is  based  upon 
the  principle  of  id  rises;  potential  energy.  She  second  variation  is  then 
taken  and  the  buckling  pressure  detersined  (Ref  1i£?5/-  the  buckling 


ores  sure  is  that  pressure  which  will  rasa! 


beinr  positive  for  a  statue  state  of  eouili 


the  secern  varramon 


The  potential  energy  of  the  system  is  given  by 

v*  =  +  ”k  *  -Q  /, 

lyvgrg  If  is  fhfc  2^2*hT32S8  3t?ci33  U.  IS  tc2  hwgir*?  st?si2 

S1A  O 

energy,  and  ft  is  the  potential  energy  of  the  external  forces.  Fro 

px  +  «I  +  -Lt^r*sldxde  , 


~  l  fr 

ih _  i  {U  4-^- 

,  _p,  i  i  f  ‘l  ^  */ 

4{i-^")i  j  L  *x  -e 

}  J 


t:b?6 


21/71  7}  +2{l -1/S /?  ”  dxde 
/X  Is  x  '  o<s 


I  L 

Li 


;  - 

- .  3w 

,,2  3u 

dwYI 

rw4 

-t-  21/ - 

4-5/  -w  - 

-r  ti  — -  ft 

L 

T  3© 

3x 

|j 

/*  —  J  >  «  L»*  4_  *  ?  I 

€x  ~  -x  T  Kx  h  2  •  x 


fl  4.  J_  l  l  4-  Ur  } 


£  =.  V  —  W  +  £r_  —  j 

^-q  &  ©  * 


7  =  Ug  +  V  +K  +  t  +fi  (1  +  p.g) 

xe 


7}  = 


LL  4-  LI 
i  ©e  i 


w  +■  £L 

X©  ;  X© 


(3-5) 

(-36) 

(3?) 

(5B) 

(40) 


?ssre 


^  *  1  '*»* 


sad  are  strains  is  the  shells  sia-pilasa  aad  %  's^s 


?}  are  the  earra.  tares  of  the  shell1  s  sLd^olsre,  la  addition,  u,vrw 
/~s 

represent  the  dissssionaHess  cosenssis  of  the  displacement  Teeter 
d*f  1  m  eg  the  eqallibriss  condition  before  hackling  occurs .  Also  /fH  C 
are  tee  -rirtaal  displacements  of  the  shall,  aquations  (55)  to  (60)  are 
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2  »  |^(27rf  *  2  2 

lKx  +He-jxf  -t  ZVKx{oa- fJL)  + 

A  L 


2  f  2  j  ,2 

+  ir?px  +Ke4'H  +  2VM/W'4 


(I  -V)^xe  +  ux  £  +  V  t  +/Aej 


+  (v6-w)  P^eJ  +v/i^  +/lx{l-v)(l,+JXe)(Lb4-vs<) 

(_  J 


+  tL  Eh'^K  +  2^  +  t2-»yt-hK)Xx  pj>axdo 


where; 


20  -  'A 

Eh  dT 


Assuming  a  pressure  distribution 


p  =  pe  +■  p,  cose 


The  results  of  a  giuphical  solution  for  equation  (69)  are  presented  in 
Figure  4.  These  curves  are  for  peak  pressures,  at  8  =  0,  which  result 
in  buckling  of  the  shell.  Thus 

P  =  P0  +  P,  (72) 

and  the  non-uniform  pressure  distribution  described  by  the  parameter  ’jfr 
defined  as  (Ref  10:1) 


3x 

w 

a 

*  m 

--  s<3 

1 
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Ejynaric  Collapse  Theory 

This  theory  is  developed  using  an  energy  analysis  and  an  approxi¬ 
mate  deflection  shape.  The  derivation  makes  use  of  the  principle  of 
minima  potential  energy  which  is  based  on  the  notion  of  virtual  work. 
This  approach  yields  a  system  of  nonlinear,  coupled  differential  equa¬ 
tions  in  the  time  varying  coefficients  of  the  deflection  shape.  Insta¬ 
bility  is  predicted  when  these  coefficients  increase  in  time  without 
bound  (Hef  3-6),  In  this  paper  instability  is  assumed  to  occur  whenever 
the  maxima  moment  applied  to  the  shell  by  the  blast  exceeds  the  critical 
maximum  moment  predicated  by  Brazier* s  analysis  (Ref  6:109). 

It  is  assumed  that  the  moment-curvature  relation  fc-r  pure  bending 
is  valid  at  each  cross-section  and  is  given  by  (Ref  3ill) 


where 


where 


M  =  -ET(7?-£?7?5) 

(74) 

g  =  L5a4(i  -VZ)j  h2 

(75) 

TL  2 

TD  =  Y  —  —  Y  Y 

~  XX  2  XX  X 

(76) 

CO 

Y(x,t)  y^An(t  ){1  coy  t  x) 
n=i 

( 77 ) 

is  the  assumed  deflection  shape  of  the  center  line  of  the  shell. 

The  basis  of  this  analysis,  as  mentioned  previously,  is  the  princi¬ 
ple  of  virtual  work  which  can  be  formulated  ass 


If  an  elastic  body  in  a  dynamic  state  is  subjected  to 
a  variation  of  Its  deformation  at  a  fixed  time,  the 
change  in  strain  energy  equals  the  incrs»3nt  of  work 
of  the  surface,  body,  and  inertial  forces  (Ref  3*H)- 


Assuring  that  there  are  no  body  forces,  the  strain  energy  can  be 
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written  as 

q/L  ,  R  4 

eg2 - ~Tj  )  dx  (78) 

'0 

Substitution  of  equation  (?6)  into  equation  (?8)  and  takiiig  the  first 
variation,  equation  (73)  becomes 


The  expression  for  the  virtual  work  of  the  inertia  forces  is 


where  u  represents  the  deformation  stats,  the  dots  represent  differ¬ 
entiation  with  respect  to  time,  is  the  mass  density  par  unit  length 
ana 


??=  =5{i  -UZ)/hZ  (91) 

Equations  (79)  and  (80)  include  terms  thru  only  the  fourth  order  in  Y 
and  its  derivatives. 

Since  only  ovslisation  terns  due  to  longitudinal  bending  have  been 
considered  and  since  equation  (73)  has  no  contribution  due  to  the  ova- 
lizatien  effect  of  the  external  pressure,  it  is  consistent  to  omit  the 
corresponding  contribution  in  the  virtual  work  expression  (Ref  3:13). 
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Only  the  virtual  work  of  the  center  line  deflection  due  to  the  external 
pressure  3s  included.  Thus 


f*  u  d3 


f 


p  co  s*e  §Yds 


ys 


where 


P  =  Prcore  e_r+  j  -  T  “  @  "  "T 


P  =0 


,  g  <33r 

/  o  “  o 


(823 


(633 


fev) 


p  is  the  peak  reflected  preaaire  am  P  is  the  decay  parameter  and 


defined  as  p^/l*  .  Substituting  (83)  into  (82) 


P  cos- 9 S Yds  =  —  Pr*e 


,-rt 


SY<1) 


>b 


(85) 


Using  a  single  tern  t^prexisation,  the  deflection  shape  is 

Z'TTx 


Y(x,+  )  =  A,(t }  (j  -  COS'  ^  } 


(86) 


Then,  by  perforating  the  appropriate  differentiation  on  equation  (86)  as 
required  in  equation  (79)  the  variation  in  the  strain  energy  becomes 


Su  =  EI 


1 6-tr  M  61 r 


74 

u 


f  -  f  *■  -fi  *,) 


!  l8  a,  I 

2 


(8?) 


the  virtual  ■gork  bee csss 

„8_2 


A  +  "bOTT  7? 
A«  l8 


§W;  = 

sa4  the  pressure  nork  becomes 


2-..  U  /  3  L 


A,A}  +  A,  A,  ]  f 

U1 


£8/ 


8w  =  |'7rpr.ae 


-rt 


^Sa. 


2  i 


using  the  psdisd.pl  8  of  TstniTwi  potential  energy 


gu  -  Sw.}  4-  gwp=  0 


(89) 
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12  8.8  E  IT4  f  C?{~}2 
+  - 1  — 

(t>4  nV  L 


tt2a,2  6 


0.375(1  — 


Equation  (93.)  is  solved  using  a  Hunge-Kuita  isethod  coupled  with  an  Mass- 
Sashforth  aM  Maas-Koulton  predictor-corrector  scheme  (Ref  12).  This 
scheme  has  provisions  for  varying  the  tdss  increment  and  allows  for 
variations  in  the  shell  sizes.  A  copy  of  this  program  is  contained  in 
Appendix  A  along  with  the  appropriate  information  n&eded  for  operation. 

The  maximum  assent,  at  ajy  given  tiwt  occurs  at  the  ends  of  the 
cylinder.  For  the  case  of  interest,  however,  the  hinge  collapse  lias 
fores  approximately  at  the  middle  of  the  shall,  ’33ms  the  marigsia  assent 
in  the  expressions  below  is  that  value  of  the  moment  at  the  middle  of 
the  shall.  Using  aquations  (7&)  thru  (77)  and  the  condition  8  «  0, 
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this  assent  is 


Recall 


where 


TTETA|L2(|  -£a*  L4  ) 


IT  (D4-d4) 


IT  0 


16(D/h  } 


D/h  >  ! 00 


d4  =  (D-h)4  =  04(!  -  £)4 


Then 


max 


1T3E  D2Affi  1.5  IT**  (D  h)2A2(t)  _ 
4  (D/h)(L/of  [(I  -  V2}  D2(WD)4 


To  detarstm  if  a  shall  amer  a  particular  olast  loading  will  fail 
crsot  involves  -fee  faHosiag  process: 

1.  Betersias  Brasher* s  collapse  sosant  fros 

_  0.4935  EO3  .  . 

Mc“  D  2  , - -  198) 

(i) 

2.  Use  fee  cossuier  program  (Appendix  A)  to  solve 
equation.  (9l)  for  the  Has  history  cf  A^_(t) 

aM  the  beading  eossnt,  H 

sax 

3.  Find  the  tigs  at  nhieh  the  sossnt*  H  , 

ST 

rushes  a  msartflua. 

4.  If  the  ttes  is  less  than  0.160  sec  and  K_,y  is 
greater  than  Brasier*  s  collapse  tangent.  Sign  the 
shell  will  coIIsdss  . 
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Static  Collapse  Theory 

Two  approaches  were  used  in  the  static  collapse  situation.  The 
first  approach  failed  due  to  the  difficulty  in  solving  the  governing 
stations  of  notion  (see  Appendix  3).  The  second  method  used  by 
Oreesspon  {Ref  10)#  yields  only  an  upper  bound  to  the  dynamic  iso-daaage 
curves.  This  fact  is  seen  in  Figure  5* 

The  first  approach  taken  utilises  the  f  orcc-equ ilibriua  sqnatiors 
fro~  shell  theory.  These  equations  are  given  by  Tisosheeko  (Ref  14)# 
and  are  generalized  to  include  the  end  effects  and  asyssetric  pressure 
loadings.  The  development  of  these  equations  for  shells  of  finite 
length  (short  as  well  as  long)  is  in  Appendix  B.  A  computer  progras  for 
shells  ox  i/D  ratios  >  15  is  also  given  in  Appendix  B. 

The  second  approach  talon  follows  that  presented  by  Greenspan 
(Ref  10).  The  pressure  loading  corresponds  to  the  pressure  parameter 
‘¥r  s=  0.5#  where  is  defined  as 

p 

'  .  .  (93) 

P  +  Pn 

Assuming  that  seshrana  theory  holds#  and  that  the  pressure  distri- 
butioK  can  he  expressed  as 


4-COJT-8)  ^  |Df{e) 

l) 


P  -  Pi  I  -t-  - uoj.  o  j 

the  gorbrane  force  resultants  can  then  be  written  as 


Ne  =apf(8] 

*  ’  o'  ■  ' 

_ „  x  af(e) 

nx6  ~  ■  °  as 


NX=  a  2 


2 

ox  -  •  A  +  Vdp_f(s)  _  -2.  -Z- 


ae2 


&  L2  a2f(©) 


*=*  24  ss2 


(\oo) 

{\  o  1) 

002.) 

(103) 
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At  the  center  of  the  shall,  :t  =  L/2, 0=0  ,  equations  (100)  thru 
(103)  .reduce  to 


P  “  Pa  +  R 


Ng  =  ap 

N  =  A 
:te  ^ 


(404 

(iOS 


L2  * 

nx=  PC-31'- T5-V)  i 

dsing  Yon  Msss  yield  condition 

:4x  -  NxH0  +  N0  +  ^*e  =  (C^h)*  i 

stisre  vJl,  Is  the  yield  stress  in  pars  tension  .\ad  h  is  the  shell 
thickness,  the  critical  pressure  for  collapse  is 


(p)c  =0^—  — — -  1  ^  - -~2  — - 

^\/(!— 2/4-1'2}  +l?/(i^2V)(-“j)  -hip  (109} 

Yor  =  0*5  sad  ^  =  0,3,  equation  (10>)  reduces  to 


(Wc  _ 


V  0/T9  -t  0.0416 (“p^)  t  0.000l085\”p_] 
the  results  of  equation  (11G)  are  shmei  graphically  is  figure  6  for 
enricas  L/2  sad  D/h  ratios. 
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V.  Psterais&tloB  of  Critical  SasH  Parsssisrg 


To  determine  f-fte  critical  shell  parcsstsrs  it  is  first  ncbesrtry  to 
construct  a  fsssly  of  iso-dansge  curves,  The  Rations  ih*t  are  us* '  to 
construct  the  iso-^iassge  are  su^arissd  first,  Use  re  rolis  of  a  atsqy 
of  these  curves  for  different-  ssteilsls  are  then  discussed. 

The  characteristic  damaged  shape  for  ths  buckling  side  is  isHcatod 
in  figure  7,  Let  the  asymptotes  of  the  tangeah-asdaliis  branch  sad  elss- 
tic  branch  be  represented,  by  !«*.  ?«  and  lg,  ?g  respectively. 

53ie  value  of  ?7»  fros  nnssricai  integration  of  equation  (16) 

(Eef  2s252)  is 

D  =  -5-  <T  I  0*0 

T  4  y  <3  r 

*h«"  o-3  is  ae  yi«M  strss£* 

Sis  value  cf  1^  can  be  estimated  by  assuring  ini  E*.  are  con¬ 
stant  during  the  plastic  flss  phase  of  bac£Ling.  Use  asilitads  of  ths 
nth*  flesnraX  sods  can  be  expressed  by  use  of  equaticr,  (24)  idier* 

”  O  .  Urns 


w(T) 
n'  * 
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n' 


[0} 


^in'nq,,T 


(u2) 


The  so  si  unpliiicd  isode  is  ths  one  having  Sis  largest  qn  aid  is  dster- 
sbK  hy  substituting  equation  (26)  into  (25)-  Thus 


Oy 


£.  \z 


'-M 

sii 


-s  max 


2E*a 


l  E 


u\t\ 

»i  * 


lihsra  is  considered  aasii  compared  to  unity  sna  has  been 

replaced  cy  o~,,  in  accordance  Kith  the  asssnstisn  cre  ranaiss  constant 
— » 

(Bef  2s23S)»  Asm *dsg  that  an  amplification  factor  of  1-000  produces 


can  be  eespressed  as  (Ref  2:236} 


f%  *tj  / *iv%  /y*  n  •.  M 
VSCL**/  mjf 


buckling,  c|nJ  max 


Bh  01 1 


C  I 


(m\ 


imdx  ”“T 

*&ere  o  is  a  constant  and  c  is  defined  as  \fi ”//>  1  .  Thas  equating 
equations  (13,3)  sad  (114)  yields 


I  =  2hBaVCTu/> 
T  3 


°y 


(US) 


%■  relaxing  the  requirement  that  remain  constant*  equation  (115) 
can  b©  brought  into  closer  agreement  with  iie  actual  situation.  Thus 
1st  S|.  be  represented  by 


<T 


Et  = 


U\Q 


K(£-6y) 

rfsre  ids  slope  of  the  stress-strain  plot  of  cr  /e  vs  percent  strain, 

£  *  defines  the  value  of  K.  By  equating  the  kinetic  energy  imparted 
by  the  impulse,  1^,  to  the  strain  energy  absorbed  in  plastic  work, 
equation  (115)  be  replaced  by 

,1  J_ 

a  (ur) 


(ff(-yfKf) 


where  B  is  assigned  the  value  12  to  match  the  numeric;  1  integration 
with  an  amplification  of  1,000  (Ref  2s238). 

Pg  is  detersdnei  from  equation  (45)  as  the  smallest  critical  value 
and  for  J/'-  3  is  (Ref  2:233)  given  by 


“  E  |x)  (x) 2 


m) 


ibe  value  of  Xg  from  (Ref  2:233)  is 
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ao  construct  the  iso-damage  curves,  the  approximate  hyperbolic 


formula. 


Wi¬ 
ll  *A 


020} 


Is  usea  (Ref  2.233)  •  Depending  on  the  curve  branch.  p  arp  j  sre 

A  A 

the  values  given  by  equations  (111)  or  (118)  and  (117)  or  (119)  respec¬ 
tively. 

(/sing  equations  ^118)  and  (119)  the  values  of  the  assgrptoies  ware 
calculated  for  1040  steel,  6061-T6  aluminum,  5052-K18  aluminum  and  1100-0 
aluminum.  These  values  are  tabulated  in  Appendix  C  for  various  L/33  and 
D/h  ratios.  Solving  equation  (120)  for  I,  then  permits  one  to  pick  a 
pressure  and  solve  for  the  corresponding  impulse.  By  picking  different 
values  for  the  pressure  and  calculating  the  corresponding  ius>ulses,  it 
is  possible  to  construct  a  iso-damage  curve  for  any  sized  shell.  Pro¬ 
ceeding  in  this  manner,  the  theoretical  curves  in  figures  12  to  16  were 
constructed,  Superiiroosed  on  these  same  Figures  are  the  experimental 
iso-damage  curvts  that  were  obtained  by  Schuman  (Ref  13:43-6), 
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la  tiis  cases  filers  Schuman  did  not  test  a  shell  of  particular  inter¬ 
est,  the  curves  os  pages  55*58  of  (Ref  13)  were  used  to  calculate  a 
critical  pressure,  The  empirical  fonaila  for  use  of  these  curves  is 
given  by 


P 

cr 


p  r  F  F.  F 
O  w  D  t  L/q 


02.*} 


where  PQ  i3  the  ataso spheric  pressure  and  the  other  functions  are  as 
described  on  Figures  8-11,  The  quantities  Fw,  F^,  F^,  and  Fj^d  are 
nondimensionai  and  are  scalars.  Using  this  value  of  critical  pressure, 
the  curves  on  pages  43-46  (kef  13)  were  extrapolated  to  give  the 
desired  ioo-daaavr  curve. 

JlgZL-es  12  to  16  present  sets  of  iso-damage  curves  for  the  four 
ssaterials  investigated  in  this  paper.  The  theoretical  curves  were 
constructed,  using  the  dynamic  buckling  theory  developed  in  chapter  4 
and  the  expression  for  the  critical  pressures  and  impulses  presented  by 
equation  (120)  in  chapter  5*  The  expreiisental  curves  ware  constructed 
using  the  ezpredmental  data  generated  by  Schucan  (Ref  13).  Below  is  a 
brief  discussion  of  each  of  the  above  figures  and  the  trend  illustrated 
by  each  figura. 
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Figure  12  is  s  set  of  iso-danage  carves  for  1040  steel.  Expert- 
mental  curve  #1,  is  for  shells  33-37  (Ref  13*97)  vfaioh  failed  in  the 
hinge  collapse  pattern.  The  experimental  p.ints  scattered  about  curve 
#2  are  for  shells  8-16  (Ref  13*97)  vhach  failed  is  the  hinge  collapse 
pattern.  These  ieo-dasage  curves  sees  to  indicate  that,  for  L/D  ratios 
of  3  snd  D/fa  ratios  in  the  range  indicated,  the  dynasic  buckling  theory 
can  predict  failure  due  to  the  formation  of  the  Msgs  collapse  pattern. 

figure  13  is  a  set  of  iso-dasage  curves  for  606I-T6  alaslnan. 

Curve  #1  is  for  shall  16?  {fief  13*100)  vMJi  failed  in  the  lings  collapse 
pattern.  Curve  #2  is  for  shell  137  vhich  failed  in  the  hinge  cell  suss 
pattern  ■with  slight  buckling.  Thus  this  shell  (#137)  failed  in  Kiat 
can  be  called  a  threshold  coronation  pat-  srn.  These  iso-damage  carves 
also  sees  to  indicate  that,  for  L/D  ratios  of  1  and  D/h  ratios  as  indi¬ 
cated,  the  dynamic  buckling  theory  again  s eess  to  predict  failure  due  to 
the  f creation  of  both  hinge  collapse  and  cgJbination  patterns. 

Figure  ltis  a  set  of  iso-da«age  curves  for  1100-0  alumnus. 

Curve  #1  is  for  shells  194-197  {Ref  13:101)  vhlch  failed  in  a  hinge 
collapse  pattern.  Curve  #2  is  for  shells  174-179  {Ref  13*100-1)  vhieh 
failed  in  a  save  buckling  pattern.  Curves  #1,  seas  to  indicate  that  for 
L/D  rating  of  3  and  D/h  ratios  of  250,  the  dynamo  buckling  theory  does 
not  accurately  predict  hinge  collapse  failure.  This  theory  scold,  how¬ 
ever,  yield  a  set  of  critical  pressures  cad  impulses  ufaicb  «col&  result 
in  an  ecicnsLve  degree  of  damage  due  to  the  forsation  of  the  hinge  col¬ 
lapse  pattern.  Carves  #2  only  seen  to  verify  that  if  a  shell  fails  in 
the  u«T3  buckling  noiec  the  dynamic  buckling  theory  is  valid  in  its  pre¬ 
diction. 

Figure  15  is  a  set  of  iso-damage  curves  for  5G52-H33  aluslrasu 
Curve  #1  is  for  shells  89-91  (Ref  13*93)  -which  fails!  in  the  hinge 
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cailapse  patters.  Curve  #2  is  for  shells  $2-$4  (Sef  13*S9)  'sbich  also 
failed  in  the  hinge  c-ollapss  pattern.  Curve  #3  la  for  shells  ^-84 
(Ref  13*93)  which  failed  is  the  combination  pattern.  Curves  #1  and  2 
seen  to  indicate  that*  for  L/B  ratios  of  2  and  3  and  for  2/h  ratios  of 
375,  the  dynamic  buekltng  theory  can  predict  failure  due  to  the  form¬ 
ation  of  the  hinge  collapse  pattern.  Curves  #3  sj»  B/h  ratios  of  500# 
the  dynamic  buckling  theory  cannot  accurately  predict  failure  due  to  the 
formation  of  the  cosbiaaticn  pattern. 

Figure  16  is  a  set.  of  iso-dasage  curves  for  5052-H3S  alasLnua. 

Curve  #1  is  for  shells  123-123  {Ref  13*99)  and  L/B  of  3  which  fsdlsd  in 
the  combination  pattern.  Curve  #2  is  for  shell  63  {Ref  13*99)  and  L/D 
cf  5  which  also  failed  in  the  combination  patters.  For  these  shells, 
note  that  they  virtually  fora  the  sans  iso-dasage  curve  and  that  they 
failed  in  the  sane  node,  Ihe  dyaagio  tackling  theory*  however,  differs 
radically  from  the  eaqperiasntal  results  and  in  addition*  differs  for  the 
L/D  ratios.  Hence,  the  dynamic  buckling  theory  cannot  accurately  pre¬ 
dict  failure  for  the  L/D  and  B/h  ratios  indicated  when  the  shell  fails 
due  to  the  formation  of  a  pattern  different  than  that  of  wave  hackling. 

Based  on  the  above  areiysi3  of  Figures  12  to  16,  it  is  concluded 
that  the  L/D  ratio  is  not  the  sost  important  parameter  in  determining 
the  failur  pattern  of  the  shells  investigated.  In  addition,  the  dynamic 
buckling  theory  can  predict  failure  due  to  the  formation  of  the  hinge 
collapse  and  combination  patterns  only  for  lew  and  moderate  B/h  ratios, 
but  can  accurately  predict  failure  for  the  wave  buckling  pattern. 

Proceeding  on  this  basis*  Figures  1?  to  23  were  constructed  by  holding 
the  L/p  ratio  constant  and  allowing  the  B/h  ratio  to  vary.  In  this 
manner,  a  limiting  value  on  the  B/h  ratios  can  be  indicated  such  that 
the  valid  range  for  the  theories  presented  in  chapter  k  can  be  established. 


this  analysis  will  also  predict  which  failure  pattern  is  to  be  sa&eoied 
given  m  Z/v  ratio*  a  B/h  ratio,  and  a  critical  reflected  pressure. 

Assuring  that  the  D/h  ratio  is  the  critical  shell  parameter  in 
d&heraihisg  the  failure  sod©,  estimates  'sere  ssde  by  this  researcher  as 
to  uhich  failure  pattern  would  appear  and  at  which  critical  pressure, 
Ihsss  shells  were  then  fired  by  Schusan  at  SHI  laboratories  is  iferyls 2d. 
Ths  results  of  these  firings  are  presented  in  Appendix  D.  For  clarity, 
however,  a  brief  discussion  will  be  presented  below, 

fhs  first  three  figures  of  Appendix  B  are  for  1100-0  Al,  The 
failure  pattern1  s  obtained  ware  as  predicted,  ie  the  collapse  pattern 
for  moderate  B/h  ratios  (275)  and  the  combination  pattern  for  slightly 
higher  D/h  ratios  (375-^0).  Ibis  information  coupled  with  data  fros 
Reference  13,  pages  10G~101,  shells  #l?b-lSi,  illustrate  the  collets 
range  of  failure  nodes.  Shells  #l?b-181  having  failed  is  the  wave  buck¬ 
ing  pattern  with  shall  parameters  of  l/D  =  1  and  B/h  *  500. 

Urs  last  two  figures  of  Appendix  D  are  for  5052-H38  Al.  Sis  fail¬ 
ure  patterns  obtained  ware  as  predicted,  is  the  wave  buckling  pattern 
for  high  B/h  ratios  (750-1000}.  Kobe  that  even  though  the  L/D  ratios 
were  radically  different  (0*333  sad  3*0),  the  sane  failure  pattern  was 
observed  for  these  high  D/h  ratio1  s. 
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As  mentioned  on  page  29 


.  a  program  was  written  to  numerically  inte¬ 
grate  equation  (91)*  The  results  of  this  program  were  not  conclusive. 

The  solution  for  A,  (t)  and  H  was  found  for  a  PESO  values  bet- 
i  max 

ween  10  and  80  ard  a  decay  parameter  (  p  )  between  9*5  and  11.1.  The 

valu,.  of  A-.  (fc)  and  M  _  thus  obtained  could  not  be  duolioated  for- 

max 

identical  input  data.  This  discrepancy  is  believed  due  to  the  inte¬ 
gration  step  size.  Hence  the  step  size  was  reduced,  but  no  output  could 
be  obtained  fro®  the  program.  By  sedifing  the  program,  interred icate 
results  were  obtained.  The  results  indicated  that  the  maximum  moment 


generated  never  exceeded  Brazier's  collapse  moment.  The  coefficient, 

JU  (t),  was  also  found  to  be  oscillatory.  Increments  on  PESO  still 
did  not  produce  a  large  enough  moment.  A  possible  explanation  for  this 
is  that  the  value  of  P  may  be  to  large.  This  large  valve  of  P  was 
assumed  such  that  regions  of  high  pressure  and  low  repulses  were  investi¬ 
gated.  The  reason  being  that  this  area  corresponded  to  a  similar  area 
generated  by  the  dynamic  buckling  theosy.  For  these  reasons,  no  theo¬ 
retical  dynamic  collapse  curves  are  presented  on  Figures  17-28 .  Perhaps 
a  value  of  P  between  0.1  and  2  0  will  yield  useful  dynamic  collapse 
curves. 

Below  is  a  short  discussion  oi  Figures  17-28  and  the  importance  of 
each.  The  curves  were  constructed  using  the  static  and  dynamic  theories 
presented  in  chapter  four  of  this  paper.  The  conclusions  drawn  from 
these  figures  are  based  on  the  assumption  that  the  limiting  value  of  the 
D/h  ratio  is  that  value  where  the  two  static  curves  and  two  dynamic 
curves  intersect  each  other.  That  is,  one  value  is  obtained  w*  ere  the 
static  curves  intersect  and  likewise  for  the  dynamic  curves. 

Figures  17-19*  which  are  for  1040  steel,  indicate  that  the  lower 
transition  value  at  B/’h  is  between  90  and  130. 
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figure  18  indicates  that  the  noser  transition  value  of 


r%  /i  » 

jjn  2.5 


approximately  435- 


Figures  20-22,  which  are 
lower  transition  value  of  D/h 


for  (06l~T6  aluminum,  indicate  that  the 
is  between  90  and  130 « 


Figures  23-25,  which  are  for  5252-H38  aluminum,  indicate  that  the 


lower  transition  value  of  D/h  is  between  90  and  135- 

Figure  24  5  dicates  tnat  the  upper  transition  value  of  D/h  is 
approximately  455. 

Figures  23-25,  which  are  for  5052-H38  aluminum,  indicate  that  the 
limiting  D/h  ratio  is  between  $0  and  130. 

Figure  2?  indicates  that  the  upper  transition  value  of  D/h  is 
approximately  960. 

Since  the  lower  transition  value  of  D/h  is  the  same  for  each  of  the 
four  materials  investigated,  assume  an  average  value  of  110.  Since  L/B 
did  not  sees  to  affect  the  lower  transition  point,  assume  that  it  does 
not  affect  the  upper  transition  point.  Hence  from  Fig-ores  18,  24,  and 
2?  this  upper  transition  value  is  approximately  at  the  value  of  D/h  = 
450.  Sow  the  failure  modes  can  be  broken  down  into  three  areas.  From 
D/h  =  0  to  110,  the  hinge  collapse  pattern  should  be  dominant.  From 
D/h  =  100  to  450,  the  combination  pattern  should  be  dominant.  From  D/fc 
above  450,  the  save  buckling  pattern  should  be  dominant. 
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VI.  Conclusions 


iihe  conclusions  reached  in  this  report  are: 

1.  Hie  D/h  ratio  is  the  dominant  shell  parameter  in 
determining  the  failure  sods  of  side-on  air  blast 
loaded  cylinders# 

2.  For  vjh  ratios  belos  300,  the  dynamic  bidding 
theory  can  approximately  predict  hinge  collapse 
and  cosoination  failures. 

3.  For  D/h  ratios  of  0  to  110*  -hinge  collapse  pat¬ 
terns  should  predominate,  110  to  453  the  com¬ 
bination  pattern  should  predominate,  and  450  up 
the  vrave  budding  pattern  should  predosdn&te. 

4.  Solution  to  the  dynamic  hinge  collapse  equations 
cannot  be  obtained  using  large  values  of  the 
pressure  decay  parameter. 
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VII,  Recosgjgmdatioss 


1*  It  is  reccmsnded  that  numerical  integration  of 
equation  (91)  be  attested  using  decay  parssater 
(  P  )  values  between  0,1  and  2,5, 

2,  it  is  recossended  that  a  series  of  ezperiisenis  hs 
conducted  to  prove  or  disprove  the  Uniting  value 
of  S/k  =  450.  Inis  sight  be  done  with  the  know¬ 
ledge  that  the  dynamic  collapse  curves  probably 
lie  below  the  static  collapse  curves  of  figures 
17-28, 

3.  It  is  recossanded  that  the  solution  to  the  equ¬ 
ations  in  Appendix  B  be  a  tiers  ted.  Their  solu¬ 
tion  would  yield  a  least  upper  bound  to  the  binge 
collapse  case. 
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Appendix  A 


Ihe  input  data  necessary  for  operation  of  this  program  is: 

1.  l/D,  D/h,  shell  diameter  (DIA)  according  to  format 
statement  100 

2.  PI  and  Poisson's  ratio'’  (NU)  according  to  format 
statement  105 

3.  pressure  parameter  (FRED),  luass  density  per  unit 
volume  (RHO),  and  Young's  modulen  (EfO)  according 
to  format  statement  110* 


The  output  of  this  program  is:  j 

| 

1.  Brazier's  moment  (MDB),  B/h,  and  L/D  ratios  j 

according  to  foraat  e-tateaent  U5  I 

| 

2*.  time  (x),  iwnriwiw  moment  (V KAX),  Brazier's  I 

moment  (MSB) ,  and  pressure  decay  parameter  (ALPHA) 

according  to  format  statement  120  I 

3.  coefficient  A,  (t)  (Y(2)#  pressure  parameter  *  1 

(FRED)  j  lounges  modulus  (ETO) ,  and  maximum 

moia?nt  (HMAX)  according  to  format  statement  125.  f 


. . mi . . . . . . . « . . . . 
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LD 
IB 
0EA 
PI 
HU 
EO 
RHO 
PRED 
d  to 
1(2) 
MSB 
X 

ALPHA 


PROGRAM  VARIABLES 


L/D  ratio 
B/h  ratio 
shell  diameter 
Pi,  3.14159 

PoissOn* s  ratio 
Young *s  modulus 
mass  density  per  unit  length 
PRESSURE  PARAMETER,  py  x  106  /  E 
Cll  coefficients  of  equation  (91) 

Coefficient,  A^(t) 

Brazier1  s  collapse  monant,  equation  (9?) 
time,  t 

dec^y  parameter,  /  1^ 


71 


non  o  n  n 


GAH/xC/68-13 


DYNAMIC  ANALYSIS  PROGRAM 

S0MM0N/CMN1/M,N»H»X,Y(2I,F(  2)»OLK2)  »PHI  (2)  ,0L0Y(2) 
REAL  Ln,LD2,Ln4,LD6,LD8, NU,MC3,MMAX 

PROGRAM  PROPER 


DIMENSION  B(l) 

READ  (5,1051  PI,  NU 

READ  (5* IOC  5  LD,DfA,ALPHA,DH,H,XM4X 
READ  (5,1105  PRED*  RHU,  EYO 
WRITE  (6 , 1C5  5  PI,NU 

WRITE  (6,100)  LOtOIA, ALPHA, DH,H,XM&X 

DO  12  l  =  l,?,l 

PJ.2  =  p  l*p  I 

p 14  =  p?2*PI2 

P  16  =  P 14*  PI  2 

PIS  =  C(4*PI4 

0H2  =  OH»DH 

DH4  =  DH?*GF2 

0 IA2  =  01 A*OI A 

ZUL  =  1.  -  NU*NU 

ZUL  2  =  7UL*ZUL 

L02  =  LC*LD 

L04  =  L02*LD2 

LD6  *  L04*L02 

LD3  a  L04+L04 

Cl  =  32»*LG8 

C2  =  C1*01A2 

C3  =  5.*DH4*ZUL2*PI 3 

C4  =  { 64.4*01 A*OH) /l.CE+12 

C6  =  (  I28.8*Pt4)/(L08*OlA2*1.0G*-06) 

C8  =  1.15*ZUL*PI4*QH2 
C9  =  L04*DIA2 

CIO  =  10.15625*®! 8*0H4*ZUL2)/L08 
Cll  =  (&.C*PI2*LO2)/0f A2 
AB  =  EYO/RHO 
C7  =  C&*AB 

MCB  =  ( (0.4935)/(0H?.*SQRT{  ZULU  ) 

WRITE  (6,115)  MCS,OK,LD,EYO 
M  =  0 
N  =  2 
Y(l)  =0. 

Y<2)  =  0. 

X  =  0. 

C5  =  C4*-A8*PRE0 
CALL  RUNGA 
IF  If)  1,2,1 

(ci/^C2  *C3*(Y(I)**2) ))*(  (C5*EX»(8ETA)  5  *  <C7*(C8*(Y( 1 )**2 
,}  I  C9)*Y(1)1  -  (C10*Y( l)*(Y( 25**2) )  +  C7*C11*(Y( i)**3) ) 

~E  (  1  >  =  Y(  2  ) 


PI,  NU 

LD,DIA, ALPHA XMAX 
PRED,  RHU,  EYO 
PI  ,NU 

LO,OIA, ALPHA, DH,H, XMAX 
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GO  TO  3 

2  8(1  )  =  V{  l } 

y ( 1)  =  Y ( 11/OIA 

MMAX  =  UPI*°I  25/(LO**2)  5*((*25/OH>-(  (  <  e  375*P  I  4*0H*ZUL) /(LD**4)  } 

1*  ( Y  {  I  )**?} J } C  1  } 

IT  (ABSlMMAX) .GT.MC8)  GO  TO  26 
Y \ 1 )  =  Btl» 

IF  (X-X«AX)  2,25,25 

25  PREO  =  »R  EC  *■  10. 

WRITE  16,125)  Y'(  l ) ,  PRE  0,  X  »  mmax  ,  ALPHA 
GO  TO  10 

26  Y  ( n  =  on  ) 

WRITE  (6,1255  Y( 1),°RED,X,MMAX, ALpHA 
OH  =  OH  +  100. 

PREO  =  I«. 

IZ  CONTINUE 

GO  TO  5 
C 

C  FORMAT  STATEMENTS 

C. 

1A0  FORMAT  (6F12.4) 

105  FORMAT  (2E15.7) 

HO  *=0RM4T  (2E15.7) 

115  FORMAT  UHC, 10X,5HMCB  =,F 1 6. 8, 5X, 4H0H  =  ,E11 .3,5X»4HL0  =*E11.3,5X, 
15  HE YO  =  1 1 1 1 *3 ) 

125  Format  (5  X,7HY(1)  =  ,E16.S,5X,6HPRED  =,E12.4,5X,6HTTHE  =,E16.8 
1,5X,7HMMAX  =  ,E12.4,5X,7HALPHA  =,E12.4) 

E-NO 

SUBROUTINE  RUNG4 

C  0M.MON/CMN1  /M,.N*H»X»Y(2)»F(2)  ,0L0  ( 2  5  ,  PHI  (2)  , OLDYf  2) 

M=M+1 

GO  TO  (15,11,12, 13,14), M 
15  RETURN 

11  DO  20  J=1,N 

20  OLDY(J)=y(J) 

DO  21  J=1,N 

21  0  K I  ( J  5  =  F  ( J  ) 

00  22  J*1,N 

22  Y(J)=OLDY( J)+.5*H*F{J) 

X=X+. 5*H 

RETURN 

12  00  23  J = 1 » N 

23  BHI { J ) =PHi { J ) +2.*F( J  5 

00  24  J  =  1,M 

24  Y  (  3  )  =OLOY  ( J  }  +-.  5*H*F  (  J  ) 

RETURN 

13  DO  25  J  =  1 , N 

25  PHI(J)=PHI(J)+2.*F( J) 

00  26  J=1,N 

26  Y(J5=0L0Y(J)+H*F{J) 

X=X+.5*H 

RETURN 
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14  no  2?  J=1*N 

27 

00  23  J*l.N 

28  YlJl=OLDYt  Jl+PH!C.J>*H/6, 
N=0 

RETURN 

END 


UAS/fS/OB-O/? 

Ippesdis  B 

Si®  governing  differential  equations,  iralraHng  *>*««-»»  ^gag©. 

ef -curvature  terns  are  given  tgr  Tisoshenko  (Ref  4:15)  as 

af?x  +  ~  <aN*«Xrv  -  NJv^r  W  .)  -  O 


c 

,  3NX  . 

f^fex  _ 

aft,'/  w 

3x 

3© 

°  X©  XX 

3N  _ 
-t-a  — x© 

+  ^©x 

_  ±£E© 

as 

3© 

a© 

aide 

N  4- 

i*sx  , 

32M^e 

*  ©  * 

3x3© 

3x2 

3©  3x 

f-dN 

xWxx  + 

2NX©(VX 

«iere 

n,  v,  w 

are  shown 

in  figure  1 

(6-3) 


are  shorn  in  figure  B-i.  lie  relation  bstsreen  ths  forces  g*d  d-igpiaf*^- 
Jssnts  and  noiasnij  and  displacsaants  are 


..  _  r 

*  ”  ~*lux 


- 


r^*[y* +  t  (v0-w)] 

£  h  r  i  ,  i 

~w^  * 

-  e  h  r  ?  i 

2(l-t-v)  a  Us  +  Ux 


L  =  ah  ...r± 

Z(\*-V)  |  a 


Mx  ~ 


WhtTC 


~Diw«  +  j  K  *  weo)j 

-o[  (V»  *Wsa)  *  V'.vJ 

-M«=  D6-^)[^(vx  +  wxj] 

>-  -^-T- 

M  (i  - V Z) 


Cb-h) 

(B-S) 

(B~c) 

( B-7 ) 

(B-S) 

(&-9) 
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Define  non-dissnsional  coordinates  as 


V 


W 


\j. 

a 

V 

a 

w 


a 

5  a 

equations  (B-l)  to  (S-3)  can  be  written  a* 

%  +  +  F- 


(B-U) 

$-n) 

(B-H) 


i-£iy-g 

L  rh 


£h  V$S 


+  Th  (¥  -wf5 


=  o 


(B-15) 


J  2~^  v||  *  Ve»  “  %  +  «*  W4|e  +  wee*  +  (l-^) 


+  V- 


*  0-^) 


Jk 

Eh 


^y|  +  ve  -w  -  3 
+  0-V1) 


'll  +  2W|S*C  * 


v*s  +  T?(Vwi\i=0 

(B-l^ 

(i~*/)-'ttB  *  Vw«  ! 


where 


§-wf|+  Tf(1  +  WH 


d=  _LfAs: 

q  \z'a) 


=  o 


£  h 


(B-17) 


(B-IS) 


*aie  ssthod  of  solution  of  these  equations  inrolTes  the  assumption 
that  the  shell  is  inextonsicnai  in  the  circnsferential  direction  or 

v0  -  w  *=■  o  (8'5<?) 

and  the  displacement  ccs§»nents  car.  be  css. pressed  as 
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M  ~  Ks  C$5  ©  4-  (|\  c«%  2$  +  *  *  . 

V  -  g.(§)  Ji.18  f  q4(|)  sin  2©  v  .  .  4 

«  -  -°'se  4  2:^(f)cn  18  4-*  ♦ 


(B-EOi 

(s~aO 

(S~2Z] 


Apoljing  equation  (3-1$)  to  the  systes  of  equations  (S-l)  thru 

(3-3)  results  in  a  set  of  equations  that  are  insensitive  to  the  drcurw 

ferential  esciensional  node,  In  the  resulting  equations,  equation  (S-l) 

resains  the  sarse,  the  linear  tera  in  K.  is  eliminated  in  the  second 

9 

and  third  equations,  and  the  linear  tern  in  S_  is  eliminated  fros  the 
resulting  expression  by  use  of  equation  (B-l),  These  manipulated  set  of 
equations  rsy  thus  be  written  as 

2  jjr  +  -  aN„Vx.  _  Me<W„„  -W)  -  O  (8-231 


a  — 4g-  -v  — —  '-L  AM©  _  _L  qZN6 
3xBe  Bx«©  3  a  3e* 


■  »» 
3x*  S© 


a  oN, 


iaN«v« 


Hr  ‘^xe  ^*e  ~^x) 


+  *-[a2 +  aw  (V  - *0l  +  alii 


?  Tn»(vo+w«) 


-*r  3UX\NX%  4-  2.  Mxe(V  -V  NA/XX){  —  d-^-^  =  Q 

J  SO1 


(&-2H) 


Assuring  P(9)  can  be  expressed  as 


p(6)  at  p  4.  p^  cos  8  ^  p  cost©  p. 


and  using  equations  (B-20)  to  (B-2X),  the  governing  set  of  equations 


become 
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hectic 


g,'  (|}  51.' 9  -  |j  C6J  5  ■%-  q,  {, 


'(f)  eo*  6 


“  sm2  si  ~  ®  (B-2&) 


,  ///t,  fi  -*  l  #  h  i c 

9i  2  cos'"  9  -*-  cos  29  —  ~  \ir)  S* 


\  ,  H  \2.  ////  _ 

—  i  —  i  «  cos  9 


h,""^}005^  + 


Si^‘~  © 


%"($)  -  vet) 

_  ii  ~;£/Z)  £  cos  s 

E(-i-)  (B-2?) 


h?w  cos  2e  +  (i-z/J 


q2'C|)  “  2h2'l^)j  Cos  2© 


—  0—  Z'-of^zCf)  —  h»f 5lt^  ZQ  ®  C*'2^ 


“2  Kj" ^  ‘  )  8  cos2  29  +  cos  49  ~  "2"  ^tnw(|)  cos  26 


k  !*H  *«  • 

-  V  1$) 


I  cos  t© 


r , 


Sin*  26 


_  pz  gps  26 

E  l"a"  )  \8-29) 

odlufioji  © t  equations  (B»26)  (B*29/  is  sot  *t~gr^~~  aejysssrrH  gng~| 

aasljrticslly,  Since  there  is  no  standard  procedure  to  solvg  thee-s 


n 


Stir/HC/&J3 


coupled  uos»l±n$ar  differential  equations  ir?  sassricsl  isiergratioa# 
th«lr  solution  is  net  sttenpied  is  this  ps>-er. 

that  if  the  sana  si^lifing  assumptions  used  is  (Rsf  4*19} 
sre  ussi  on  equation  Cis-2o)  thru  (B-29),  the  resulting  expressions  are 
the  sag®  as  found  os  page  i$  of  (Hef  4). 
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Appssiis  G 


Sie  input  data  necessary  for  aeration  of  this  pregras  is: 

1.  Young* s  sodnitts  (SO)  and  3*S3  dsnsiiy  (REG) 
aacordiTig  to  format  statassgt.t  100 

2.  sotting  the  initial  value  of  ShriL  radius,  A 

3.  salting  the  initial  value  of  i/j  ratio*  IS 
4„  setting  the  initial  3/h  ratio 


Sie  output  of  this  prograa  is: 

1,  Young*  s  ssaslus  (SO)  and  sass  csnsiiy  pa?  unit 
vdisss  (HHO)  according  to  fsrsat  statcssat  135 

2,  the  asymptotic  Tgluss  of  pressure  and  ispuLse 
(rS  aid  IS),  rail ns  to  length  ratio  {AL}* 
diaster  to  i.hicsness  ratio  (HA),  and  pressure 
parameter  (fH25)  according  to  fornat  stateast 

115 

3,  the  pressure  (?)  srcl  irpuiso  (IK?)  for  iso- 
-dasage  carves  of  the  elastic  scdcdus  hackiii^ 
theory  according  to  the  fcrsat  ststosit  120 
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PROGRAM  VARIABLES 


p 

pressure*  ?r 

PE 

asEypototic  pressure 

HIED 

pressure  parsneter  PE  x  10 fy  E 

IMP 

iepuise,  Ij. 

IE 

asymptotic  iapulse 

no 

Young*  s  nsodulus,  E 

RHO 

isass  density  per  unit  volume 

A 

shell  radius 

LD 

L/D  ratio 

AL 

radius  to  length  ratio 

HA 

a/h  ratio 

-SI  8  JOB _  MAP  — 

tIP.FTC  MAIN 

.  _ REAL  IE  _t_,LD 

REAL  IMP 

'  5  READ .  < 5 , 10_0±_EYO  *  RHO 
WRITE  (6,105)  EYOi RHO 

.  _ A  _ 

LD  =  X. 

_  BETA,j=  (RHO*EY05/J12t*32.2) 

DO  11  J  =  1,3,1 

...  At  r  4LP  .  . 

HA  =  50. 

_ P.D..10.  '  xl„ _ 

PE  =  (0.92*EY0*AL)/(HA**2.5) 

„.  ..  ...  IE.  =.  ( 5.*A*SQRT  ( BETAJ  *1000. ’  / (HA**25 
•  PRED  =  (PE*1.0E+6)/£YO 

WRITE  (6,115)  AL,HA*PE_»!E, PRED 
P  =  PE  +  5. 

_ DO  12  I  =  l_»10_»l _ 

IMP  =  (  (  TE*PE) /(P-PEl)  +  IE 
___  WR ITE.  (  6 , 1 2 0J_  P_,J MP 

p  =  p  +  io.' 

12 _ CONTINUE,  _ 

HA  =  HA  +50. 

10  CONTINUE _ 

LD  =  LD  +  2. 

11  CONTINUE 
GO  TO  5 

100  FORMAT  (£16.8,J12.4) 

105  FORMAT-  (lHd’»r15X»5HEYO  =,E16.S,lbX,'5HRHO'  =,El2*4)  '  ^ 

115  FORMAT  (1H0»5X«4HAL  ^»E10.2,5X,4HHA  =»E12.4,5X,4hPE  =,E16.5,5X, 
1  AH IE  =,£16.8  *5X, 6HPRED  =,£16.8) 

...220  ..  FORMAT!  1HC,3HP  =,E16.8,5X_»5KIMP  =,E16.85 

END  .  . . 
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TABLE  C  -  6 
MATERIAL  PROPERTIES 


1  / 
1  ,a>Sk'V'<. 

E 

ibf/m2 

P 

ifam/ih2 

£ 

K 

!bf/m2 

|  1040 

6 

7< 

|  STEEL 

3  Ox  10 

.258 

34,500 

6061-T6 

ALUMINUM 

10.0  x  10° 

.098 

30 

40,000 

5052-H  38 

5 

ALUMINUM 

10,2  x  10 

.097 

32 

37,000 

1100-0 

6 

ALUMINUM 

10.0  x  10 

.098 

34 

5,000 

*  Kef  (  2 :  21?) 
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Appendix  D 


The  following  table  of  blast  parameters  is  for  the  shells  appearing 
in  this  section. 


TABLE  0  -  1 
BLAST  PAKAESTERS 


FlCSJRE  0'3 

Coshinsties.  1100  -  A1 
0  isi..  L/D  =  3.4.  D/h 


fh  *  400 


TtKinHhUlUUilHHlUliKtltf  jMillliliill 


mk 


Richard  Gonrad  Fees  ws s  bom  on  H  August  i 945  in 
Pittsburg*  Pennsylvania.  He  graduated  froa  Franklin  Area  High 
School  in  Marxysville,  Peraisylvarda*  attended  the  University  of 
Pittsburgh5 s  Greensburg  campus  until  Jamary  1965  and  then 
transferred  to  the  sain  cassus  in  Pittsburgh.  While  at  the 
min  caffijus,  he  was  elected  to  Pi  Tan  Sigm  National  Mechani¬ 
cal  Ssginsering  Fraternity,  Sigma  Tan  National  Engineering 
Fraternity,  Outstanding  Sophsore  in  Engineering,  and  Oo~Chair~ 
aan  of  Engineer1  s  Week  1?6?*  He  also  was  a  nszsber  of  the 
Arnold  Air  Society  and  Scabbard  and  Elsie  Military  Honorary 
Society.  He  graduated  in  August  1$6?  Ctua  Lands  and  received 
Ms  Bachelor  of  Science  in  Mechanical  Engineering  and  a  cos. 
mission  in  the  USAF.  His  first  assigssent  was  residence  at 
AFIT  -  SE. 

Permanent  address:  4146  03d  Ek  Penn  Hay. 

Hurrysville,  Pa.  15668 

Ibis  thesis  was  typed  by  Mrs.  Mary  E,  Pees, 


"IlllitllU1 


«I.  SUPPLEMENTARY  NOTE* 


i  to  special  export  controls  sad  each 
nationals  say  be  sade  only  with 
U  Air  Force  Institute  of  Technology  (AFIT-SJ 
&33 _ " 

U.  tPOHUHUHt  Ml  .lTAflY  ACTIVITY 


13.  ABSTRACT 

~  The  purpose  of  this  investigation  is  to  determine  the  influence  of  the 
size  paraseters,  length- to-dianater  ratio  and  diaseter- to- thickness  ratio*  on  the 
failure  sods  of  cylinders  under  side-on  air  blast  loading.  These  failure  nodes 
are:  il)  hinge  collapse,  £)  wave  buckling,  andU3)  co^dnation.  results  of  this 
study  show  that  the  disaster- to- thickness ,  (D/h),  ratio  drnrfTiates  the  gh»n 
behavior  and  that  for  values  of  D/h  below  110,  hinge  collapse  -soass  to  be  dosinants 
between  210  and  4^50 ,  combination  failure  seess  to  be  dominant;  above  45®*  wevs 
buckling  scans  to  be  doidnaat.  l 
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